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ABSTRACT 

The high mass X-ray binary 4U 1901+03 was reported to have the pulse profile 
evolving with the X-ray luminosity and energy during its outburst in February- 
July 2003: the pulse peak changed from double to single along with the decreasing 
luminosity. We have carried out a detailed analysis on the contemporary phase- 
resolved energy spectrum of 4U 1901-1-03 as observed by Rossi X-ray Timing 
Explorer {RXTE). We find that, the spectra are phase dependent. At the begin 
of the outburst, the maximum of the optical depth for Compton scattering is 
near the major phase peak. During the decay of the outburst, the optical depth 
has the maximum being away from the main peak of the pulse profile. For each 
observation, Fe Ka emission line is detected in the phase-resolved spectra, and 
its flux is constant across the pulse phases. This suggests an origin of Fe emission 
from the accretion disk but not the surface of the neutron star. 

Subject headings: pulsars: individual (4U 1901+03)-stars: neutron-X-rays: stars 

1. Introduction 

Most of the known X-ray binary pulsars are the so-called high mass X-ray binary 
(HMXB). They are usually the Be/X-ray binaries characterized by the transient nature. 
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The Be star is an early-type nonsupergiant star with observable emission lines from the ma- 
terial in its circumstellar disk (see Slettebak 1988 for a review). Be/X-ray binaries usually 
show up with two types of outburst behavior: normal outburst with low X-ray luminosity 
lasting for days- weeks, and giant outbursts with higher X-ray luminosity (Lx ^ 10^^ erg 
s~^) occurring irregularly in every several years. Giant outbursts are thought to be driven 
by a dramatic expansion of the disk surrounding the Be star, which leads to the formation 
of an accretion disk around the compact object. Accompanying with the giant outburst, 
the neutron star is detected with the pulsed emission. The spectra of Be/X-ray binary are 
usually represented by a cutoff power law shape (e.g., Coburn et al. 2002; Corbet et al. 2009; 
Crawford et al. 2009). The iron K feature between 6 and 7 keV and low-energy absorption 
due to the cool material are as well observable (White et al. 1983; Wilson et al. 2008). 

4U 1901-1-03 was detected as a Be/X-ray binary pulsar, one member of the HMXB (Liu 
et al. 2006). The orbital period and eccentricity of the system are measured as 22.58 d and 
0.035, respectively. The period of X-ray pulsation is about 2.73 s (Galloway et al. 2005). 
The location of the source is R.A.=19''04™13^4 and Dec=+3°09'26" (J2000.0) obtained by 
the observation with Uhuru (Forman et al. 1978; Priedhorsky & Terrell 1984). With the 
observation of i?XTi? Proportional Counter Array (PCA), Galloway et al. (2003a) obtained 
precise coordinates of R.A. = 19''03™37M, Dec=+3''ll'3l". 

There are two giant outbursts in the history of observations of 4U 1901+03. The first 
one was detected with Uhuru and Vela 5B in 1970-1971 (Forman et al. 1976; Priedhorsky & 
Terrell 1984). The second giant outburst took place in 2003 February, and was first detected 
by All Sky Monitor (ASM) on RXTE, followed by series of pointed RXTE observations over 
the next five months (Galloway et al. 2003b). At the hard X-rays, contemporary observations 
from INTEGRAL satellite are also available. During the second giant outburst, the X-ray 
fiux of the source reached a value of -F2.5-25fcey~8xl0~^ ergs cm~^ s~^. 

Based on the RXTE observations. Galloway et al. (2005) made a thorough study on 
the orbital parameters, preliminary X-ray spectral analysis and pulse profiles from several 
observations for 4U 1901-1-03 during the 2003 giant outburst. Chen et al. (2008) carried 
out the detailed analysis of the energy dependence of the pulse profile along the outburst, 
and found that the pulse profile is correlated with both the X-ray luminosity and photon 
energy (Wang & Welter 1981; White et al. 1983; Nagase 1989; Mukerjee et al. 2000). 
The phase-resolved spectral analysis is need for studying the emission configure of Be/X-ray 
binary pulsar. Here we present for the first time phase-resolved spectra of 4U 1901+03 using 
all available RXTE data taken during the 2003 giant outburst. The paper is organized as 
the follows: description of observations and spectral model in Sec. 2, results in Sect. 3 and 
discussion finally in Sec.4. 
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2. Observation and data reduction 

The observations analyzed in this paper are from PCA and High Energy X-Ray Timing 
Experiment (HEXTE) on board the RXTE satellite from February 10 to July 16 2003. The 
PCA consists of 5 non-imaging, coaligned Xe multiwire proportional counter units (PCUs) 
covering a nominal energy rang from 2 to 60 keV. Only PCUO and PCU2 data are adopted 
in this work, which were on all the time during the observation of 4U 1901+03. The HEXTE 
instrument consists of two independent clusters (cluster A and cluster B) covering an energy 
range from 15 to 250 keV. Because detector 2 of cluster B of HEXTE lost its spectral 
capability and automatic gain control, only cluster A data are analyzed in the work. We 
extract light curves and spectra of PCA and HEXTE from intervals when the source has the 
offset angle of less than 0.02° and the limb of the earth is more than 10° with respect to the 
source direction. All HEXTE data products are dead-time corrected using the HEASOFT 
ftool hxtdead. 

The data of Standard-2 and GoodXenon modes of the PCA and Archive and Science 
Event modes of the HEXTE are used to perform the spectra and timing analysis. The total 
lightcurve is extracted with Standard-2 mode data, and pulse-averaged X-ray spectra are 
extracted with Standard-2 and Science Event mode data. The pulse profiles of the chosen 
observations are extracted with GoodXenon, and the phase-resolved spectra are extracted 
with GoodXenon and Science Event data modes, using the software f asebin. Figured] shows 
the lightcurve of 2.0-21.0 keV, the soft and hard colors that are defined as the count rate 
ratios 4.5-6.1 keV/2.0-4.5 keV and 9.8-21.0 keV/6. 1-9.8 keV, respectively. 

The PCA background subtraction is carried out using the latest versions of the appro- 
priate background models, and a 1 % systematic error is added to the spectra to account for 
the calibration uncertainties. Events in energy range ~ 2.5-20 keV of the PCA and 17-80 
keV of the HEXTE (17-50 keV for the late observations of the outburst) are selected for 
the spectral analysis with the software XSPEC version 12. 3. Op (Arnaud 1996; Dorman & 
Arnaud 2001). 

2.1. Spectral Model 

According to the characteristic of the lightcurve (Fig. 1), we study the phase-averaged 
spectra of ten typical observations, carried out on February 10, 15, 22 and 23, March 27 and 
30, April 30, May 25, June 14 and 28, respectively, , covering the outburst profile of the 
beginning, the peaking, the stepping down at the middle of decaying and the ending tail. 
Various spectral models are used to fit the phase-averaged spectra of these observations. 



-4- 



e.g., cutoffpl, powerlaw, bknpower, compTT, and a combination of each with a blackbody 
component. A gaussian component centered at 6-7 keV which represents fluorescent Fe hne 
emission shows up as well in all spectra. A model consisting of compTT with a spherical 
geometry (Titarchuk 1994) and a Gaussian component is statistically acceptable {Xred ~ 1) 
to flt the spectra (also see Galloway et al. 2005). Due to low effective area PGA below 3 
keV, the column density nn of neutral absorption can not be constrained well. Therefore, 
for all spectral fits we fix nn at 1.2 x 10^^ cm~^ (Galloway et al. 2005). Such a model can 
well fit the data till April 30. For the observations near the end of the outburst (May 25, 
June 14 and 28), an additional component of a blackbody has to be introduced in order to 
have a reasonable x^ed- "^^^ temperature of the blackbody is A;Tbb ~ 1 keV. As an example. 
Figure [2] shows the fitting results with various models to the data of April 30. 

We choose three typical observations on February 10, March 27, and April 30 to analyze 
the evolvement of the phase-resolved spectra. The phase-resolved spectra are extracted with 
eight phase bins. Figure [3] shows the spectra from the data of April 30, for phase bin No. 5 
(0.4375-0.5625). The fitting results of the spectra show the gaussian component is needed 
in the spectrum. 

3. Results 

3.1. Variation of soft and hard colors 

Figure [1] shows a stepping down feature of the light curve, happened at time around 
April 30, during the decay of the outburst. The diagnosis on the pulse profile revealed a 
change from double peaks to single peak at this time (Ghen et al. 2008). Accompanying 
to this as well are modifications in the fit model of the spectrum (Galloway et al. 2005; 
Ghen et al. 2008). Along the decay of the outburst, the hard and the soft colors have the 
different evolutions: the soft color evolves similar to the light curve, while the hard color 
shows the opposite and ends up with a rapid decrease. The soft color and the hard color 
suggest an overall trend that the spectrum softens at the begin, but hardens with the decay 
of the outburst, and turns to soft at the end. 

3.2. Phase-averaged spectra 

Table 1 shows the results of the fit models on ten typical observations described in 
Sect. 2, i.e., the beginning, the peaking, near stepping down and near the ending. The 
results show the temperature of the seed photons (To) in the Gomptonization model is about 
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1 keV, and decreases slightly during the decay of the outburst. The temperature of the hot 
electronic population (kT) remained at about 5 keV. The optical depths of the observations 
from February 10 to April 30 are consistent with a flat distribution within the error bars, 
and an average is derived as ~ 5. At the end of the outburst, the optical depths are obtained 
as 7-8. We notice that, the Fe Ka line emission with equal width ~ 100 eV presents at 
around 6.5 keV in the energy spectrum, and its flux varies between ~1-10 x 10^^^ ergs cm~^ 
s~^, showing a trend of increasing with the intensity. At the end of outburst, an additional 
spectral component of a blackbody is needed for spectral fittings. 

3.3. Phase-resolved spectra 

For studying the evolution of the phase-resolved spectra, we choose the typical observa- 
tions with the high signal-to-noise ratio, i.e., on February 10, March 27 and April 30. The 
spin lightcurve of each observation is subdivided into 8 phases, for each the energy spectrum 
is extracted, and the flux is shown in Figure HI Figure H] shows the phase dependence of 
the spectra of the three observations. The pulse profile of the three observations are best 
presented in Figure 4, each with 32 phase bins. At the beginning of the outburst (February 
10), the optical depth maximizes near the major phase peak. However, during the decay 
of the outburst, such maximum moves towards the second one, and the temperature of the 
scattering electron has a trend of anti-correlation with the optical depth. The center energy 
of Fe line is almost around 6.5 keV. For each observation, the fluxes of iron line are likely 
constant across the whole pulse phases. The spectra from the phases (1 for February 10, 2 
for March 27 and 2 for April 30) of the lowest intensity present as well the existence of the 
obvious Fe feature, as shown in Figure is an example from the phase 2 of the data of April 
30, with the fluxes comparable to those from what else phases (Fig. H]). 

3.4. Spectral Ratios 

The detailed phase-to-phase variation is best illustrated by ratios of the phase-resolved 
spectra to the spectrum of the phase with minimum count rate (Leahy & Matsuoka 1990). 
We show the phase-resolved spectral ratios of the three observations. The minimum-count- 
rate spectra are from phase 1 for observation of February 10, phase 2 for both observations 
of March 27 and April 30. The phase ratios are then presented separately in Figure. [MHl 
There are dips in the PHA plots at ~ 6.5 keV (as marked by the arrow in the first panel 
of Fig. ED. The dip is due to that flux of Fe emission line that only occurs in non-pulse 
component, this is consistent with the results of the phase- resolved spectra. 
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For the phase 2 of February 10, at the beginning of the main pulse, the ratios increase 
with energy. The ratio increases with energies rapidly at the phase 3 where the main pulse 
is peaking, but slowly at the phases (4 and 5) beyond. As is the case as well for at the end 
of the main pulse (phase 6), till the energies around 10 keV, beyond which the ratio remains 
almost constant. Such a trend holds more or less in phases 7 and 8, corresponding to the 
second pulse, but with a turn over of the ratio at energies larger than 10 keV. For March 
27, the properties of PHA ratios are similar to those for February 10. For April 30, phases 
3, 4 and 5 correspond to the main pulse peak. The PHA ratios of these phases are increase 
with energy, but are constant at above 10 keV. Phases 6, 7, 8 and 1 are around the second 
pulse, and their PHA ratios increase at low energy (<10 keV), and decrease at 10-20 keV. 
The overall spectral ratio of April 30 is different from those of February 10 and March 27. 

The PHA ratios of the phase-resolved to the phase with minimum count rate are in- 
creasing with energy, illuminating that the spectra of the pulsed are generally harder. The 
stronger the flux of the pulsed, the harder the spectrum. At the phases around the second 
pulse, the spectrum decreases abruptly at higher energies. 

4. Discussion 

The X-ray radiation modes of Be/X-ray binary pulsar is generally thought to have 
tight relation with the luminosity (Parmar et al. 1989). At high luminosity (> 10^^ ergs 
s""^), the accretion flow onto the magnetic pole can be decelerated via the radiative shock 
formed near the neutron star surface (Wang & Frank 1981). The emitting plasma will be 
compressed under the shocked region, where the photons can only escape from the sides of the 
column (fan-beam mode). Under some circumstances, a pencil-beam may still emerge (Nagcl 
1981). At lower luminosity (< 10^'' ergs s~^) the infalling material may be decelerated in a 
coUisionless shock above the neutron star surface (Basko & Sunyaev 1975; Kirk & Galloway 
1981). In such a case, a thin emitting region or effects of the strong magnetic field can cause 
a pencil-beam of emission to be formed (e.g., Meszaros et al. 1983). The various models 
could be verified by the observational properties of X-ray pulsars. 

We have analyzed the phase-resolved spectra of 4U 1901-1-03, and found that the op- 
tical depth and the temperature of the scattering electron are related to the pulse phases, 
which is a common feature of X-ray pulse binary (e.g.. La Barbera et al. 2003). Spectral 
ratios indicate the main pulse peak has the hardest spectrum, which is common property of 
accreting pulsars (Hickox & Vrtilek 2005; Tsygankov et al. 2007). Our results show, at the 
beginning of the outburst (February 10) where the luminosity (> 10^^ ergs s~^), the emission 
of the main pulse has the possible origins of the fan-beam. The main pulse has the large 
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optical depth, could be due to that, during the main peak, the angle between the column axis 
and the observer's line of sight has the highest value so that the observer is looking almost 
along the beam (Klochkov et al. 2008). During the decay of the outburst, our results are 
consistent with the emission configuration that the main peak from the fan-beam and the 
second peak from pencil-beam (Chen et al. 2008), the main peak that has the hard spectrum 
is due to that high energy photons are more liable to escape in a fan-beam from a hot region 
close to the footstep and perpendicular to the accretion column (Basko & Sunyaev 1976; 
White et al. 1983; Klochkov et al. 2008), which is consistent with the low optical depth and 
high temperature of the scattering electron of the main peak. The second peak is produced 
by the low energy photons escaping in a pencil-beam along the direction of the accretion 
column where the optical depth is lager and the electron temperature is lower. These results 
are accordant to those of Chen et al. (2008) who analyzed energy-resolved pulse profiles in 
details. They concluded as well that, the fan-beam contributes to the main pulse peak and 
pencil-beam to the second pulse peak. In addition, in some energy bands the fiux varies by 
almost a factor of two between different pulse phases (Chen et al. 2008), both the angles of 
the spin-axis from the viewing direction and between the magnetic pole and the spin-axis 
are substantial. Therefore, it is likely that emission from both magnetic poles contribute to 
the pulse profile. 

We notice that the evolution of hard color is similar to that of pulse fraction (see Chen 
et al. 2008) , and spectral ratios show that the spectra of the pulsed are harder than those 
of the non-pulsed under larger pulsed fluxes. This indicates that the high energy photons 
contribute mostly to the pulsed emission, which is consistent with the evolvement of the 
pulse fraction shown by Figure 4 of Chen et al. (2008). Perhaps the lower energy radiation 
is dominated by photons from the accretion column walls that hit the neutron star surface 
and are reprocessed, which would lead to a lower pulse fraction. 

Fe Ka emission line with equivalent widths of several hundred eV is usually detected in 
X-ray pulsar, may be caused by illumination of neutral or partially ionized material in the 
accretion disk, stellar wind of the high-mass companion, or material in the line of sight or in 
the accretion column (Pravdo et al. 1977; Basko 1980; Nagase 1985; Paul et al. 2002; Naik 
et al. 2005). Our results show that prominent Fe emission with equivalent width ~ 100 eV 
is detected in the spectra of 4U 1901-1-03. The flux of Fe emission fine is increasing with 
the X-ray luminosity along the outburst (Fig. 1 and Table 1), and is proportional to that of 
continuum intensity and increases with increasing of the accretion rate (Suchy et al. 2008). 
For each observation, Fe line is detected in the spectra of the continuum in each pulse phase 
and its fiux is constant within errors, and is not correlated with pulse phases. Therefore, the 
Fe line of 4U 1901-1-03 is not likely to origin from the accretion column of magnetic polar cap. 
Furthermore, the spectral fit of the phase-resolved spectrum shows that, the center energy 
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of Fe neither change among the pulse phases nor evolve along the outburst of 4U 1901+03. 
We therefore believe that, during the outburst of 4U 1901+03, Fe emissions should come 
from a region in the accretion disk but not the surface of the neutron star. 
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Table 1: Fit pEirameters for the spectra. 
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Note. — To is the temperature of the soft seed photons for the Comptonization, kT is the electron 
temperatiire. t is the optical depth of the scattering cloud in a spherical geometry. E'fo is the central energy 
of Gaussian emission line, and Fluxpe is the flux of Gaussian emission line in units of 10~^^ ergs cm~^ s~^. 
fcTbb is the temperature of blackbody, and norruhh is the normalization of blackbody. 
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Fig. 1. — From top to bottom: the total light curve, the soft color, and the hard color, 
respectively. 
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Fig. 2. — PCA spectrum for the observation of April 30 and the best-fit model. The 
residuals are shown sequentially taking (a) the cutoffpl model (x^g(j=3.96), (b) the pow- 
erlaw+bbody model (x^gj=13.15), (c) the powerlaw+bbody+gaussian (Xred— 6.38), (d) 
the bknpower model (x?ed=9-04), (e) the bknpower+bbody model (x^<,(j=2.09), (f) the 
bknpower+bbody+gaussian (Xrcd=l-^5), (g) the compTT model (Xred— 1-69)) (h) the 
compTT+gaussian model (x^^^—O-QS). 
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Fig. 3. — The fitting results of tfie pliase-resolved spectrum of the phase 5 of the observation 
April 30, the left panel shows the fitting result by the model without Gaussian component, 
the residual shows obvious Fe line profile; the right panel shows the fitting result with the 
model containing the Gaussian component, the residual is removed. 
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Fig. 4. — The phase-resolved spectral parameters as a function of pulse phases. Prom top to 
bottom: the pulse profiles of 32 phase bins obtained from all PCA channels (2-60 keV), the 
flux of 3-30 keV of 8 phase bins in units of 10~^ ergs cm""^ s~^, optical depth for Compton 
scattering, the temperature of the scattering electron cloud, the central energy and flux of 
Fe emission line. The central energies of Fe line are flxed for phases 4 and 5 of March 27 
and phases 3 and 4 of April 30. 
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Fig. 5. — The fitting result of the spectrum of the phase 2 of the observation April 30. The 
left panel shows the fitting result and residual, the right panel shows the model components 
of (compTT+gaus) . 
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Fig. 6. — Ratios of intensity in narrow energy bands from February 10 phase 2-8 spectra to 
intensity from phase 1 spectrum. 



-17- 



1.2 

1.15 
1.1 
1.05 
1 

1.3 
1.2 

O 1.1 

I- 

< 

DC 1.6 
< 

Q. 

1.4 

1.2 
1.4 
1.3 
1.2 
1.1 



1^ 



phasel 
phase2 

I I I I I I 



phases 
phase2 



+ ++++ 
+++++ 



j + 1 I I I I I I I I I I I 



phase4 



phase2 +^+H-+ ++ |tt± 

+ .4- + 



~ -1++*+++ 

-t^ I I I I I I I I I I I 




phases +f+/tt+,l|%- 



phase2 ^+_^+ ++I 



ft "I 



+++ 



10 15 
Energy (keV) 



10 15 
Energy (keV) 



Fig. 7. — Similar to Figure 6, but for March 27. 



-18- 



phasei 

- phase2 


' ' 1 




J 1 1 1 1 1 1 1 

\ phases 

phase2 i 1 

i-^-i-i 1 1 1 1 1 1 


(/t|,t| 




- phase4 hi 
" phase2 1 if '|l 








1 1 1 1 1 1 1 1 1 1 1 1 

- phases ,i | i, , - 

- phase2^ /||||'tt|| 

- +/++++++ ^ 
"++ 

J. , 1 , , , , 1 , , , , 1 


/I If 



5 10 15 

Energy (keV) 

Fig. 8. — Similar to Figure 6, 
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for April 30. 



